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Abstract 

We present a toy model for five-dimensional heterotic M-theory where bulk three-branes, orig- 
inating in 11 dimensions from M five-branes, are modelled as kink solutions of a bulk scalar field 
theory. It is shown that the vacua of this defect model correspond to a class of topologically distinct 
M-theory compactifications. Topology change can then be analysed by studying the time evolution 
of the defect model. In the context of a four-dimensional effective theory, we study in detail the 
simplest such process, that is the time evolution of a kink and its collision with a boundary. We find 
that the kink is generically absorbed by the boundary thereby changing the boundary charge. This 
opens up the possibility of exploring the relation between more complicated defect configurations 
and the topology of brane-world models. 
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1 Introduction 



The single most important problem in trying to make contact between string-/M-theory and low-energy 
physics is probably the large number of degenerate and topologically distinct vacua of the theory. It 
is usually stated that non-perturbative effects will eventually lift most of this degeneracy. However, 
despite the advances over recent years in understanding non-perturbative string- and M-theory there is 
very little indication of progress in this direction. In fact, with the advent of M-theory and concepts such 
as branes and brane-world theories new classes of vacua have been constructed and, as a consequence, 
the degeneracy problem has perhaps grown even more serious. It seems worthwhile, therefore, to ask 
whether the cosmological evolution rather than inherent non-perturbative effects of the theory may play 
a prominent role in selecting the vacuum state. Indeed, it is known that the degeneracy of some vacua 
(particularly among those with a large number of supersymmetries) will not be lifted non-perturbatively, 
suggesting cosmology will have some role to play. 



The first task to tackle, in this context, is the formulation of a workable theory capable of describing 
a number of topologically different vacua and transitions among them. As a second step, one will have 
to analyse the cosmological evolution of this theory. It is precisely these two problems which will be 
the main topic of the present paper. 

The class of vacua we will use in our approach is provided by compactification of heterotic M- 
theory ^ on Calabi-Yau three folds 1^, ^, ^, |5| resulting in five-dimensional brane-world theories @, ^• 
These theories are defined on a space-time with two four-dimensional boundaries corresponding to the 
fixed planes of the orbifold /Z2 and, in addition, may contain bulk three-branes which originate from 
M five-branes wrapping two-cycles in the Calabi-Yau space |2|, ^. The associated effective actions are 
five-dimensional gauged = 1 supergravity theories in the bulk coupled to four-dimensional N = 1 
theories residing on the two boundaries and the three-branes. The prospects for particle-physics model 
building within this class of compactifications is quite promising and a number of models with attractive 
particle-physics properties on the "observable" boundary have been constructed The simplest 

way to characterise topologically different compactifications from the viewpoint of the five-dimensional 
effective theories is by using the charges ai and 02 on the boundaries and the three-brane charge 
03. These charges are not independent but must satisfy the cohomology constraint ai + 02 + 03 = 
which follows from anomaly cancellation. Two five-dimensional effective theories with different sets of 
charges (01,02,03) originate from topologically distinct compactifications. A transition between two 
such theories may occur through a small-instanton transition [l^, 17 1 when a three-brane collides with 
one of the boundaries. The three-brane can then be "absorbed" by the boundary and, correspondingly, 
the boundary charge is changed by the amount carried by the incoming three-brane. This change in the 
boundary charge indicates a more dramatic transition in the boundary theory. For example, the gauge 



group and the amount of chiral matter ||18|] may be altered as a consequence of the internal topology 
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change. 



The goal of this paper is to find a five-dimensional (toy) model which provides a unified description 
for the above class of topologically distinct vacua, in the simplest setting, and allows for transitions 
between them. While, for simplicity, we will assume that the topology of space-time both in the 
internal Calabi-Yau space and in the orbifold direction remains unchanged we will allow for transitions 
corresponding to a topology change in the internal gauge-field instantons on the boundaries and a change 
in the number and charges of three-branes. Our basic method will be, starting with five-dimensional 
heterotic M-theory in its simplest form, to model the three-branes as topological defects (kinks) of 
a new bulk scalar field x- We do not claim, of course, that this model provides the correct definition of 
M-theory in these backgrounds. However, we do show that the defect model in the background of its 
various vacuum states reproduces the five-dimensional M-theory effective actions with different charges 
(Qi,a2,«3)) corresponding to topologically distinct M-theory compactifications. 

Time-evolution of the defect model and the scalar x in particular then allows for a transition between 
these topologically distinct configurations. We will study in detail the simplest such transition, namely 
the collision of a three-brane kink with one of the boundaries. This will be done by calculating the 
four-dimensional effective action for the defect model in the background of such a kink. As we will 
see from this four-dimensional action, the collision process indeed generically leads to an absorption 
of the kink and a change in the boundary charge by the amount carried by the kink. Hence, we have 
established the existence of one of the elementary topology-changing processes in our defect model. This 
opens up the possibility, subject of ongoing research, that a study of more complicated configurations, 
such as brane-networks, will provide insight into topological properties of brane-world models. 

The plan of the paper is as follows. In the next section, we will introduce the five-dimensional 
effective actions from heterotic M-theory, in their simplest form. For later reference, we will also review 
the associated four-dimensional effective theories. Section 3 then presents our defect model and explains 
how, precisely, it is related to the M-theory actions. In Section 4, we will compute the four-dimensional 
effective action for the defect model in the background of a kink and Section 5 presents the resulting 
evolution equations. Section 6 is devoted to a detailed study of the kink evolution and its collision with 
a boundary, based on these equations. A conclusion and outlook is presented in Section 7. 

2 Effective actions from heterotic M-theory 

To set the scene, we will now describe the five-dimensional brane world theories for which we would 
like to find a smooth defect-model. These brane-world theories can be viewed as a minimal version of 
five-dimensional heterotic M-theory [Q. For later purposes, it will also be useful to review the four- 
dimensional effective action associated to these brane-world theories. 
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Coordinates for the five-dimensional space M5 are denoted by where a, /3, • • • =0,1, 2, 3, 5. We 
also introduce four-dimensional indices /x, z/, • • • = 0, 1, 2, 3. The coordinate y = is compactified on an 
orbi-circle /Z2 in the usual way, that is, by first compactifying y on a circle with radius p and then 
dividing by the Z2 orbifold action y —y. Taking the y-coordinate in the range y € [— vrp, np] with 
the endpoints being identified the two resulting four-dimensional fixed planes (boundaries), denoted 
by Ml and M|, are located at y = and y = vrp, respectively. Such a geometry is obtained by 
compactifying 11-dimensional heterotic M-theory on a Calabi-Yau space. If the five-branes present in 
the 11-dimensional theory are included in this compactification they lead, upon wrapping a two-cycle in 
the Calabi-Yau space, to bulk three-branes in the five-dimensional brane-world theories. For simplicity, 
we will consider a single such three-brane whose world- volume we denote by M|. We also need to include 
the Z2 mirror of this three-brane with world- volume Mf. Three-brane world- volume coordinates are 
denoted by a^. In the minimal version of the model the bulk fields consist of the metric and the dilaton 
$ while the three-brane world- volume fields are simply the embedding coordinates X" = X"(o"'^). The 
effective action for these fields is then given by 



Note, that the dilaton $ measures the size of the internal Calabi-Yau space which is, more precisely, 
given by fe*, where f is a fixed reference volume. It relates the five-dimensional Newton constant K5 
to its 11-dimensional counterpart k via 



Further, Oi, where i = 1,2,3 are the charges on the orbifold planes and the three-brane, respectively. 
They are quantised and can be written as integer multiples 




(2.1) 



(2.3) 



of the unit charge a defined by 



o = — 




(2.4) 



These charges satisfy the important cohomology condition 
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(2.5) 



1=1 
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which follows from anomaly cancellation in the 11-dimensional theory. The quantity a which appears 
in the above bulk potential is a sum of step-functions given by 



a = ai^(Mi) + a29{Mi) + as (^^(M|) + 9{Mi)j . (2.6) 
Finally, the induced metric on the three-brane world-volume is defined as the pull-back 

7^, = d^X^d.X^g^p (2.7) 

of the space-time metric. 

For positive three-brane charge, > 0, the above action can be embedded into a five-dimensional 
= 1 bulk supergravity theory coupled to four-dimensional = 1 theories on the boundaries and the 



branes. The details of this supergravity theory have been worked out in Ref. |2C]. In the case of an anti- 
three-brane, that is for 03 < 0, while bulk supersymmetry is preserved everywhere locally, it is broken 
globally. Technically, this happens because the chirality of the four-dimensional supersymmetry pre- 
served on the three-brane is opposite to the one on the orbifold fixed planes. Such non-supersymmetric 
heterotic models containing anti-branes have not been studied in much detail, so far. We have included 
this possibility here because it will naturally arise later in our discussion of the defect model. The 
generalisation to include more than one three-brane is straightforward. It simply amounts to replacing 



the Nambu-Goto type three-brane action in the third line of (2.1) by a sum over such actions (with 
generally different three-brane charges) and modifying the cohomology condition (|2.5| ) and the definition 
of a, Eq. (2.6), accordingly. 



Note that two actions of the type ( p.l| ) but with different sets of charges correspond to topologically 
different M-theory compactifications. Specifically, the charges ai and 02 on the boundaries are related 
to gravitational and gauge instanton numbers. If we keep the topology of the Calabi-Yau space fixed, as 
discussed, different values of ai and 02 indicate a different topology of the internal gauge bundles. As 
a consequence, the values of ai, are also correlated with other properties of the boundary theories, 
such as the types of gauge groups and the amount of chiral matter. Different values of 03 imply different 
internal wrapping numbers for the five-branes and, hence, clearly indicate different topologies. 



For the case of a three-brane (rather than an anti-three-brane) , the action (2.1) has a BPS domain- 
wall vacuum 20 1 given by 

ds^ = alhdx^dx''r]^u^hlh^dy^ (2.8) 

e* = 60/1^ (2.9) 

X^ = af" (2.10) 

X^ = y = const (2.11) 
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Here the function h = h{y) is defined by 

%) = -^|«^"^" + ^° for o<M<y 

3 [ (ai +03)1^1 - asT + co for 1" < |y| < vrp 

and ao, bo and cq are constants. Note that this solution is not smooth across the three-brane reflecting 
the fact that the three-brane as described by ( |2.1| ) is infinitely thin. Such a static BPS solution does 
not exist for the anti-three-brane since the sum of the tensions ai + a2 + jasl does not vanish for < 
by virtue of the cohomology condition ( |2.5[ ). In fact, solutions which couple to an anti-three-brane will, 
in general, be time-dependent. 

For later reference, it will be useful to discuss the four-dimensional effective action associated to the 



(2.13) 



brane-world model and the above BPS vacuum. It is given by |21, 20| 

The three scalar fields </), /3 and z have straightforward interpretations in terms of the underlying higher- 
dimensional theories. The field c;/), as the zero mode of the five-dimensional scalar <I>, specifies the volume 
of the internal Calabi-Yau space averaged over the orbifold. More precisely, this average volume is given 
by ve'^ . The scalar /3, on the other hand, originates from the (55)-component of the five-dimensional 
metric and measures the size vr/je^ of the orbifold. Finally, z represents the position of the three-brane 
and is normalised to be in the range z £ [0, 1] with the endpoints corresponding to the two boundaries of 
five-dimensional space-time. The four-dimensional Newton constant Kp is related to its five-dimensional 
cousin by 

4^1;. (2.14) 

Finally, the three-brane charge 

Q3 = vr/jas = eo/33 , /?3 e Z (2.15) 



is quantised in units of eo as defined in Eq. ( p.4D and is positive for the case under discussion. 

As expected, the action ( p. 13 ) can be obtained from an = 1 supergravity theory by a suitable 



truncation. The Kahler potential for this supergravity theory has been first given in Ref. |21]. An im- 
portant quantity which governs the validity of the effective action ( |2.13| ) is the strong-coupling expansion 
parameter 



eoef'-'^ . (2.16) 



It measures the relative size of string loop corrections to the four-dimensional action ( p. 13 ) or, equiv 



alently, the strength of the warping in the orbifold direction from a five-dimensional viewpoint. The 
effective action ( ^.13 ) is valid as long as e < 1 and can be expected to break down otherwise. Another 
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reason for a breakdown of the four- as well as the five-dimensional effective theory is the five-brane 
approaching one of the boundaries, that is, z — > or z — > 1. In this case, the underlying heterotic 
M-theory may undergo a small-instanton transition |16, ^] which leads to the M five-brane being con- 
verted into a gauge-field instanton (or, so called gauge five-brane ||22|) on the boundary. In such a 
process, properties of the boundary theory, such as the gauge group and the amount of chiral matter, 
can change dramatically as a result of the internal topology change [^] . In our simple five-dimensional 
model such a modification of the boundary theory is indicated by a change in the boundary charge 



ai or a2 by the amount of incoming five-brane charge. It is clear, however, that the actions (2.1) or 
( |2.13D are not capable of describing such a jump in the boundary charge in a dynamical way. In fact, 
the four-dimensional action ( |2.13| ) does not retain any memory of the presence of the boundaries as 
z 0, 1. This can also be seen from the moving-brane solutions to ( 2.13| ) found in Ref. and will 
be explained in more detail later. As we will see, our defect model, to be presented in the next section, 
will considerably improve on these points. 



3 Modelling heterotic brane-world theories 



We would now like to find a "smooth" model, replacing the five-dimensional action (2.1), where the 
three-brane is not put in "by hand" but, rather, obtained as defect solution to the theory. Such a model 
should have, as a solution, a smooth version of the BPS domain wall ( |2.8| )-( prT^ ). Note, that we will 
not attempt to find a smooth description for the orbifold fixed planes. Their nature, as part of the 
space-time geometry, is entirely different from the one of the three-branes. In particular, the fixed plane 
tensions ai, 02 can be negative whereas the three-brane tension |a3| is always positive. 

Modelling co-dimension one objects such as our three-branes is usually achieved using kink-solutions 



of scalar field theories |1£]. This is indeed what we will do here. We, therefore, supplement the bulk 
field content of the five-dimensional theory by a second scalar field x- For this bulk scalar along with 
the dilaton <I> and the five-dimensional metric, we propose the following action 



1 
2^1 



+ 



-g 2W 



-g 2W 



(3.1) 



We require that the potential V be obtained from a "superpotential" W following the general formula |24] 

(3.2) 



V = ]-G^^diWdjW - ^W^ 



where Gij is the sigma- model metric and indices I,J,... label the various scalar fields <I>^. For our 
specific action (^.Ij), we have two scalar fields (^>^) = (^,x) &iid the sigma-model metric is explicitly 
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given by 



G = diag(^,e-*) . (3.3) 



Further, we propose the fohowing form for the superpotential 

W = e-*u;(x) (3.4) 



where w is an, as yet, unspecified function of x- Using the general expression ( p. 21 ) this resuhs in a 
potential 



1 —2(1? 2 ^ — <I>^^ / dw\ 



Note that, in ( |3.lD , we have omitted the Nambu-Goto type action for the three-brane corresponding 
to the third line of the M-theory effective action ( |2.lD . The reason is, of course, that we would like to 
recover the three-brane as a kink-solution of the new scalar field x- For this to work out, the potential 
U has to have a non-trivial vacuum structure. In fact, since the original three-brane charge is an 
(arbitrary) integer multiple of a certain unit, we need an infinite number of equally spaced minima. 
More precisely, we require the potential U satisfies the following properties : 

• [/ is periodic with period v, that is U{x + v) = U (x) 

• U has minima aX x = Xn = nv for all n £ Z 

• U vanishes at the minima, that is U{xn) = 0. 

These requirements can be easily translated into conditions on the function w which determines the 
superpotential. Clearly, from the second and third condition, the derivative of w has to vanish at all 
minima Xn = nv oi U . The definition (^]^) of w in terms of U involves a sign ambiguity which allows 
one, using the first condition on U above, to make w periodic as well. However, the structure of the 
action (^) makes it clear that the "vacuum values" w{xn) of w have to reproduce the charges on the 
orbifold planes. We, therefore, define w as 

X 



w{x) = / dxVJM) (3.6) 

JO 

which implies quasi-periodicity, that is, 

w{x + v) =w{x) +w{v) . (3.7) 

We have plotted the typical form of U and w in Fig. |l[ 

For much of our discussion the concrete form of the potential will be irrelevant as long as the above 
conditions are met. A specific example, however, is provided by the sine-Gordon potential 

'2vrx^ 



U = rr? 



1 — cos 



(3.8) 
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Figure 1: Shown is the typical shape of the superpotential w and the potential U (in units of a) as a 
function of the scalar field x ('^''^ units of v). 



where m is a constant. The associated superpotential w is easily obtained by integration. 

This concludes the set-up of our model. Let us now discuss how, precisely, this model corresponds to 
the brane-world theory ( |2.1| ) introduced earlier. The simplest solution for x is to be in one of its vacuum 
states, that is, X = Xn = nv for some integer n, throughout space-time. In this case, the superpotential 
and potential reduce to 

W = e-^wixn) , V = ^e-^'^wixnf • (3.9) 



Substituting this back into the action (|3.1| ) and comparing with the M-theory result ( |2.1| ) shows that this 
precisely corresponds to a situation without a bulk three-brane. In particular, one concludes that the 
boundary charge ai has to be identified with the value w{xn) = w{nv) = nw{v) of the superpotential 
at the respective minimum Q This is, of course, the more precise reason why we have required the 
superpotential to be quasi-periodic rather than periodic. Furthermore, we learn that the elementary 
unit of charge a in the M-theory model (see Eq. (|2.4|)) corresponds to w{v), that is, 



w{v) = / diy/uix) ■ (3.10) 







^Note that, in the absence of three-branes, we have 0:2 = — ai from the cohomology condition (2.5). Therefore, also 
the charge on the second boundary is correctly being taken care of by our model. 



The next more complicated solutions are kinks where the scalar field x interpolates between two of 
its minima as one moves along the orbifold direction. Due to the cross-couplings in the action ( |3.1| ) also 
the dilaton $ and the metric necessarily have a non-trivial profile in this case. To find such solutions, 
an appropriate Ansatz is provided by 



$ = $(y) 

X = x{y) ■ 



(3.11) 
(3.12) 
(3.13) 



The four y-dependent functions A, B, <I>, x are subject to the second order bulk equations of motion to 
be derived from the first line in ( |3.1| ) and the boundary conditions 

(3.14) 
(3.15) 
(3.16) 





= -1VF = 

3 


1 

— e w 
3 




dW 




-B 1 

e X 


dx 


dw 
dx ' 



Here, the prime denotes the derivative with respect to y and the equations hold at both boundaries, that 



is, at ?/ = and y = up. The first equality in each equation is easily derived from (3.1) including the 
boundary terms while the second one follows from inserting the explicit form of the superpotential (3.4). 

Instead of dealing with the second order equations to obtain explicit solutions it is much simpler 
to consider the first order BPS-type equations. Their existence is guaranteed by the special form of 
our scalar field potential V as being obtained from a superpotential [^]. Concretely, inserting the 
Ansatz ( p.ll[ )-( ^.13| ) into the bulk part of the action (3.1) one obtains an energy functional 



E 



4A 



dye 
dye'^ 



-6e 



-2B A/2 



1 



4 



1 



x' +v 



dW 
9¥ 



X Te 



dW 
dx 



^(3e-^^'±VF)' 



(3.17) 



y=np 
y=0 



which can be written in Bogomol'nyi perfect square form. This leads to the following first order equations 



-B 



A' 



e-^x! 



1 -* 
3 



±6" 



,dW 



dw 



(3.18) 
(3.19) 
(3.20) 



dx dx 

Again, the second equality in each line follows from inserting the explicit superpotential (3.4). The 
scale factor B is, of course, a gauge degree of freedom and can, for example, be set to a constant. It is 
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clear then that equation ( |3.20| ) for x decouples from the other two. This x equation is, in fact, exactly 
the same first order equation one would derive for a single scalar field x with potential [/ in a flat 
background. It is, therefore, clear and can be seen by direct integration, that this equation admits kink 
solutions where x interpolates between a certain minimum x = x„ = of C/ at y — > — oo and one of its 
neighbouring minima at y ^ +00. More precisely, for the choice of the upper (lower) sign in Eq. ( p. 201) 
the minimum at x = + l)f^ (x = ~ 1)^) is approached for y +00. The corresponding solutions 
for A and ^> can then be obtained by inserting this kink solution and integrating Eqs. ( |3.18 ) and (3.19). 



In the next section, this will be carried out in a more precise way. In addition, the solutions obtained in 
this way have to satisfy the boundary conditions (|3.14| )- (|3.16 ). Clearly, this is automatically the case if 



the upper sign in the first order equations ( |3.18D -( p.20| ) has been chosen, that is, if the kink interpolates 
between the minima x = and X = (^ + 1)^ foi' increasing y. For the lower sign, on the other hand, 
there is no chance to satisfy the boundary conditions and, hence, no solutions of the type considered 
here exist in this case. The interpretation of these results is straightforward. While both types of kinks 
are on the same footing as far as the bulk equations are concerned the boundary conditions distinguish 
what should then be called an anti-kink, interpolating between x = ^^t^ and X = (^ ~ 1)^) from a kink, 
interpolating between x = '^'^ and X = (^ + 1)^- While the latter represents a BPS solution of the 
theory, the former carries the wrong orientation to be compatible with the boundaries and, in fact, 
will only exist as a dynamical object. This is in direct analogy with the properties of three-branes and 



anti-three-branes in our original M-theory model ( 2.1 



For the case of a kink, we would like to make this correspondence with the M-theory model more 
precise. Let us consider a kink solution to Eqs. ( |3.18| )-( |3.20| ) and ( p.l4| )-( |3.1(]| ) with the kink width 
being small (compared to the size of the orbifold) and the core of the kink sufficiently away from the 
boundaries. In this case, the profile for x and w{x) can be approximated by a step-function. Specifically, 
we have w{x) — nw{v) to the left of the kink and w{x) — (n + l)w{v) to the right. Inserting this into 
the equations ( 3.18| ), ( |3.19| ) and the boundary conditions (|3.14| ), ( p. IS] ) for A and $ and solving the 



resulting system precisely leads to the BPS three-brane solution given by Eqs. ( |2.8D , ( |2.9D , ( 2.12| ). The 
charges appearing in this solution are given by 

ai = no" , 02 = —{n + l)a , 03 = a , (3-21) 

where we have used our earlier identification ( |3.10| ) of the superpotential value w{v) with the elementary 
charge unit a. Hence, our model allows for a solution which can be interpreted as a smooth version of 
the M-theory domain wall coupled to a single-charged three-brane. 

More generally, we would like to discuss the relation between the action (|3.1|) in the background 



of a kink solution and the M-theory action (2.1). To do this, we should allow for fluctuations of the 



kink. It is well-known [25| that, for sufficiently small width, the hypersurface prescribed by the kink's 
core is a minimal surface and is, therefore, adequately described by a Nambu-Goto action. Practically, 
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this implies that the kinetic term for x ^-^id the U potential term in the action ( [3.1| ) can be effectively 
replaced by a Nambu-Goto action describing the dynamics of the core of the kink. Of course, this core 
has to be identified with the three-brane in the M-theory model. It is easy to show that, by virtue of 



Eq. ( 3.10 ), the tension in this effective Nambu-Goto action is given by a which is the correct value for 
a single-charged three-brane with /Js = 1. Further, the superpotential w in such a kink background can 
be effectively replaced by a step-function, as discussed above. Using the identification ( 3.2l| ) of charges. 



it is easy to see that the superpotential w precisely equals the function a , defined in Eq. (2^), in this 
limit. As a consequence, the second potential term in ( |3.1[ ) proportional to e~^*tt;^ precisely reproduces 
the bulk potential in the M-theory action ( |2.lD . Similarly, the boundary potentials in ( ^.1| ) match the 
boundary potentials in ( ^ ) using that w{x{y = 0)) ~ mr = ai and w{x{y = vrp)) ~ (n + l)cr = —02- 
Although there are no BPS anti-kink solutions, it is clear that a similar argument can be made for the 
action (^]^) in the background of an anti-kink leading to the M-theory action ( |2.1| ) with an anti-three- 
brane. 



In summary, we have seen that the action ( |3.l| ) in the background of various vacuum configurations 
of the field x reproduces different versions of the M-theory effective action (p.l|). For a constant field x 
located in one of the minima of [/, we have reproduced the M-theory action without three-branes. For a 
kink (anti-kink) background with sufficiently small width away from the boundaries we have obtaining 
the M-theory action with a single-charged three-brane (anti-three-brane) . Note that, while from the 
viewpoint of the smooth model ( |3.lD these cases merely correspond to different configurations of the 
field Xi they represent different effective actions on the M-theory side. As we have discussed, these 
different effective actions arise from topologically distinct compactifications of the 11-dimensional M- 
theory. While these compactifications are known to be related by topology-changing transitions such as 
small-instanton transitions these processes cannot be described by the action (2.1). What we have seen 



is, that our smooth defect model incorporates a number of these topologically distinct configurations 
within a single theory and, may, describe transitions between them as the scalar field x evolves in time. 
In the subsequent sections, we will study the simplest example for such a transition, namely the collision 
of a kink with one of the boundaries. 

A final comment concerns the question of multi-charged branes. Clearly, multi-charged BPS three- 
branes with > 1 are allowed in the M-theory model (2.1). However, our defect model (p.l|) does 



not have exact BPS multi-kink solutions as long as the potential U is smooth at its minima. The 
reason is that, for smooth U , a kink solution does not reach a minimum within a finite distance, as 
can be easily seen from Eq. ( 3.20| ) with U expanded around a minimum. As a consequence, single- kink 



solutions cannot be "stacked" to produce exact multi-kink solutions. There are a number of options 
available to remove this apparent discrepancy. Firstly, the model ( |3.1| ) as stands does have approximate 
multi-kink solutions (with exponential accuracy) which could be identified with multi-charged three- 
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branes. Secondly, if the potential U is continuous but non-smooth at its minima a kink solution can 
reach a minimum within a finite distance. There is no obstruction then to build up exact multi-kinks by 
stacking single-kink solutions. Thirdly, some multi-scalar field models are known to admit multi-kink 



solutions [^6[. So, we may generalise the action (3T) by adding more than one scalar field. For the 



purpose of this paper, we will not implement any of these possibilities explicitly but, rather, focus on 
single-kink solutions in the following. 

4 The four-dimensional effective action of a kink solution 

We would now like to study one of the simplest dynamical processes in the context of our defect 
model, namely the time-evolution of a kink solution and its collision with a boundary. For a sufficiently 
slow evolution this can be conveniently studied in the context of the four-dimensional effective theory 



associated to (3.1) in the presence of a kink. The purpose of this section is to compute this effective 
four-dimensional theory. As we will see, this computation can be pushed a long way without specifying 
an explicit potential U . We will, therefore, keep U general throughout this section. An explicit example 
for U will be studied in the next section. 

Our first step is to write the kink solution in a form which makes the dependence on the various 
integration constant (which will be promoted to four-dimensional moduli fields later on) as explicit as 
possible. We find that the kink solution to ( |3.18| )-( ^l20| ) and ( p.l4D -( p36| ) interpolating between the 



minima X = Xn = and x = Xn+i = {n+ l)v for increasing y can be cast in the form 

X = C7(eV(y/vrp-^)), (4.1) 
= e<^(l + eoe^-'^/(y,/3,z)), (4.2) 

A = Aq + ]:^, (4.3) 
b 

5 = /3, (4.4) 

where we recall that A and B are the scale factors in the five-dimensional metric as defined in Eq. (|3.11j )- 
( p.l3D . The functions C and / in the above solution can be expressed in terms of the potential as follows 

C~\X) = ^ r , (4.5) 

/(y,/3,z) = r dywfcfe^i^-'iy/Trp-z))) . (4.6) 

Trpeo Jyo ^ ^ ' ' 

Here 0, /3, z, integration constants, while is a constant which measures the width 

of the kink in units of vrp. It is clear from the form of the metric ( |3.11| ) that the constant can 
be absorbed into the four-dimensional metric. As we will see, it is, however, convenient to keep this 
constant explicitly since it can be used to canonically normalise the four-dimensional Einstein-Hilbert 
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term. For our subsequent discussion, let us define the average {h) of a function h = h{y) over the 
orbifold by 

1 r-TTp 

{h) = — dy h{y) . (4.7) 



-^P Jo 

Since the constants t/q 4' really describe the same degree of freedom, we can fix y^ by requiring 
that (/) = 0. With this convention, the integration constant (j) has a clear geometrical interpretation, 
namely e*^ represents the orbifold average of the dilaton e*. Similarly, measures the orbifold size in 
units of Tip. The final integration constant z specifies the position of the kink's core (the position where 
X = (f^ + \)v) in the orbifold direction. Values z G [0, 1] imply that the kink's core is located within 
the boundaries of five-dimensional space and is, hence, physically present. Further, z ^ 0, 1 indicates 
collision of the kink with one of the boundaries. For z ^ [0, 1] the core is outside the physical region 
and we can merely think of z as the virtual position of the core were space-time to continue beyond 
the boundaries. In this case, the physical part of the kink, located between the boundaries, is only its 
tail. In the limiting case z — > ±00 the kink disappears completely and we approach one of the trivial 
vacuum states of the theory with either x = ^^'W or x = ('^ + throughout five-dimensional space-time 
depending on whether z — > -|-oo or 2; — > —00. Also note that the function C, defined in Eq. ([4.5| ), is 
independent of all integration constants and can be computed for a given potential U . 

We should now promote all integration constants in our kink solution to four-dimensional moduli 
fields. This leads to three scalar fields (0^) = {(f), (3, z) and the four-dimensional effective metric Qii^u- 
Accordingly, the Ansatz ( 3.11| )-( [3.13 ) should then be modified to 



ds' = ^^^y^^'^dx^^dx''g^^, + ^''^y^^'Uy\ (4.8) 

$ = $(y,,/.^), (4.9) 

X = xivA'). (4.10) 

where A, ^ and x ^re as in Eqs. (|T|)-(|;|) but with {(j)^) = (</>, /3, z) now viewed as functions of the 
external coordinates x^. 

We are now ready to compute the four-dimensional effective action. Inserting the Ansatz ([4.g|)-(4.10) 
into the action (3.1) and integrating over the orbifold direction we obtain the following result 



- / 



The sigma-model metric G/j is given by 



Gij = 2(e 



JA+B 



1. , 1 



(4.11) 



(4.12) 



-mAdjA - Zd^^iAdj^B + -di^dj^ + -e-^dixdjx 

where dj = and {(j)^) = {(l),(3,z). Further, in order to obtain an Einstein- frame action we have 
required that 

(e2^+^) = 1 . (4.13) 
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This indeed fixes the constant Aq in Eq. ( [4.3[ ) to be 



The four-dimensional Planck scale is defined by 



2ttp ' 



(4.14) 



(4.15) 



as usual. 



The remaining task is now to evaluate the expression (|4.12| ) for the moduli-space metric using the 
kink solution (4.1)-(4.4). This leads to fairly complicated results, in general. There is, however, an 
approximation suggested by the original M-theory model which simplifies matters considerably. As 
discussed, the effective actions for heterotic M-theory in Section 2 are valid only if the strong-coupling 
expansion parameter 



eoe' 



P-<t> 



(4.16) 



is smaller than one. We are, therefore, led to compute the moduli-space metric ( 4.12D in precisely 
this limit which corresponds to small warping in the orbifold direction. Concretely, we will keep terms 
up to 0{e) and neglect all terms of 0{e^) and higher in our computation. This implies a dramatic 
simplification since the function /, which enters the kink solution Eq. ( ^.21 ) with an 0{e) suppression, 
drops out at this order. Inserting (^)-( |4^ ) and ( 4.14| ) into ( 4.12| ), one then finds for the moduli-space 
metric 



G 











f + e- 



V 



e-'^((5.x)') ) 



(4.17) 



Using the solution (^ 

G/313 
Gzz 



for X we finally obtain 
1 

2' 
3 



J2(eV'(l-^)) 
Ji (eV'(l-^)) -^1 



eoe 



JQ[e^^Ji-\l-z)]-Jo[-e"^l-'z 



(4.18) 
(4.19) 
(4.20) 
(4.21) 

(4.22) 

where we recall that the function C, defined in Eq. ( fl.Sj ), can be computed for any given potential 
U and is, by itself, independent of the moduli. The above result, good to 0(e), for the sigma model 
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as the only non-vanishing components of G. Here, the functions Jn are defined by 

Jn{x) = r dx x-U{C{i)) = - dx {C'\x)T ^'(X) , 

eo Jo Jc(o) 



metric explicitly displays the complete moduli dependence of G and its only implicit features are the 
dependence on the potential U and a simple integral thereof. We find it quite remarkable that the 
calculation can be pushed this far without an explicit choice for the potential U . 



The result (4.18)~( |4.21j ) suggest the existence of another expansion parameter besides e, namely the 
quantity e~^fi. It represent the ratio of the kink's width and the size of the orbifold. Working in a 
thin-wall approximation where this ratio is much smaller than one our results simplify even further. 
Clearly, we then have to good accuracy 



3 

2 



f3z 







(4.23) 



For the remaining non-trivial component G^z we can explicitly carry out the integral ( 4.22 ) and find by 
inserting into Eq. ( [4.21 ) 



G, 



eoe^-'^F{(3,z) 



(4.24) 



where 



F{P,z) 



w 



C feVHl -^))]-w{c (-eV^^ 



- [w{y = vrp) - w{y = 0)] 



a 



(4.25) 
(4.26) 



Here, the notation w{y = 0) {w{y = irp) indicates the value of the superpotential evaluated for the kink 
solution at the boundary y = (y = irp). 

To summarise, in the limit of both the strong-coupling expansion parameter and the ratio of wall 
to orbifold size being smaller than one, that is. 



eoe 



/3-0 



<1, 4<1, 



(4.27) 



the moduli-space metric for the kink solution is well-approximated by 



G = diag( ^,|eoe^-<^F(/3,z; 



with associated four-dimensional effective action 
1 



^4 



(4.28) 



(4.29) 



Here, the function F is as defined in Eq. ( [4.25| ). 

It is interesting to compare this four-dimensional effective action to its counterpart ( p. 13 ) obtained 
in the M-theory case. Obviously, the only difference arises in the kinetic term for z where the function 
F appears in ( 4.2S| ) but not in the M-theory result ( 2.13| ). A detailed comparison requires computing 
this function from Eq. ( 4.25| ) by inserting an explicit potential U. However, the qualitative features 
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of F can be easily read off from the alternative expression ( 4.26 ). It states that F is the difference of 
the superpotential on the two boundaries in units of a and, hence, it is simply the "charge difference" 
between the two boundaries. Suppose, that the kink's core is well within the physical space and away 
from the boundaries, so that z S [0, 1] and sufficiently different from the boundary values 0, 1. The field 
X will then be very close to the minimum x = nv at the y = boundary and very close to the minimum 
X = (n + l)v at the other boundary. The charge difference between the boundaries and, hence, the 
function F, is, therefore, very close to one. If, on the other hand, the virtual position of the kink's core 
is at z > 1 (z < 0) and sufficiently away from the boundary, x '^ill be close to the minimum x = i^v 
{x = {n + l)v) on both boundaries. Hence the function F is approximately zero in this case. This 
obviously implies a non-trivial behaviour of F close to the boundaries for z ~ and z ~ 1. As a result, 
for the kink being inside the physical space and away from the boundaries by a distance large compared 
to its width the effective action ( 4.29| ) completely agrees with the M-theory result ( |2.13| ) . Conversely, 
if the kink approaches one of the boundaries or collides with it, that is, z 0,1, the function F 
becomes non-trivial and the effective theories ( [1.29 ) and ( 2.13 ) differ substantially. It is clear then, that 
the effective theory ( 4.29] ) carries some memory of the presence of the boundaries while the M-theory 
action ( p. 13 ) does not. For this reason, studying the collision process in the context of (4.29) is an 
interesting problem which we will address in Section 6. 



5 An explicit example 

In this section, we consider the explicit example of the double-well potential 

U = m\v^-X^f, (5.1) 

where m is a constant. As stands this potential does, of course, not satisfy our periodicity requirement 
for U. However, for our purposes this is largely irrelevant since the single-kink solution in which we are 
interested here probes the potential only between the two minima 0. The associated superpotential is 
given by 

w = mx (^v"^ - ^X^^ ■ (5-2) 

Hence the elementary charge unit a and eo take the form 

4 4 , 

a = w{v) — w{—v) = -mv , eo = vt/oct = -TTpmv . (5-3) 
o o 



We recall that our kink carries a single charge and we should, therefore, set /Ja = 1 in Eq. ( [2.15D to obtain perfect 
agreement. 



3 



One way to satisfy all earlier requirements is to restrict the potential (5.1) to the interval [— ViV] and continue it 



periodically outside. The subsequent results do not depend on whether one works with this periodic version of the 



potential or simply with its original form (5.1). 
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The kink-solution for this potential is of the general form ( |4.lD -( ^^ with the functions C and / 
given by 



C{x) = V tanh(x) 



(5.4) 



and 



/ 



eoe 



1 



v'^ tanh^ ^ ln(cosh ^) 



/i \7rp 



where the thickness of the kink can be identified as 

1 



(5.5) 



(5.6) 



mvTTp 

The constant c in Eq. (|5.5| ) has to be fixed so that < / >= 0, as discussed before. This leads to an 
expression involving di-logarithms and we will not carry this out explicitly. 

Instead, we consider the limit where the strong-coupling expansion parameter e remains small, so 
that / becomes irrelevant and our general result ( 4.18| )~( 4.21| ) holds. The functions J„ can now be 
explicitly computed inserting the potential (|5.1| ) and ( |5.4| ) into their definition ( 4.22| ). This leads to 



Jnix) 



dx 



cosh^ X 



(5.7) 



This, together with Eqs. ( [4.18D -( fl.21 ) completely determines the moduli-space metric for the double- 
well potential as long as e < 1. While the above integrals can be carried out for all relevant values 
n = 0, 1, 2, the cases n = 1 and n = 2 lead to somewhat complicated expressions, the latter involving a 
di-logarithm. However, Jq takes the relatively simple form 



Jnix) 



1 sinh X 

— tanh X + 



4 cosh X 



(5.8) 



As is clear from the general case discussed in the previous section, for a kink with small width, that 
is, e~^p < 1, fortunately Jq is the only relevant function. In this limit, the moduli-space metric is. 



therefore, given by the general form (4.28) which we repeat for convenience 

1 3 



G 



diag( -,|,eoe^-MA^) 



The function F, defined in Eq. ( 4.25| ), now takes the explicit form 

F{P,z) = Jo(e''fi'Hl-z)) - Jo 



(5.9) 



(5.10) 



where Jq is given in Eq. ( |5.8| ). Inserting this result into ( 4.29| ) completely determines the four- 
dimensional kink effective theory for e < 1 and e~^/i < 1. The function F above indeed has the 
properties mentioned in the previous section, namely F ~ 1 for z well inside the interval [0, 1] and 
-F — > for z — > iboo. The typical shape of -F as a function of z is shown in Fig. 
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Figure 2: The function F which enters the effective four- dimensional action of the kink as a function 
of z for e^ = 10. 

6 Kink evolution equations 

We will now study the time-evolution of the kink based on the effective four-dimensional action derived 
in the previous section. The collision of the kink with one of the boundaries will, of course, be of 
particular interest. 

We focus on simple time-dependent backgrounds and a metric of Friedmann-Robertson- Walker form 
with flat spatial sections, that is 



dsi 



(6.1) 
(6.2) 



where {(p^) = {(f>,P,z). Let us first review the general structure of the evolution equations for back- 
grounds of this form. From the general sigma-model action ( [4.11| ) one obtains the equations of motion 

^" ''^^ (6.3) 

(6.4) 
(6.5) 



3d2 



2a + 3d2 



+ 3d(j)^ + r 



0, 



where F is the Christoffel connection associated to the sigma-model metric Gjj and the dot denotes 
the derivative with respect to time. Adding the first two equations, we obtain an equation for the scale 

18 



factor a alone which can be immediately integrated. Discarding trivial integration constants one finds 

a = - In |t| . (6.6) 

This power-law evolution with power 1/3 is as expected for a universe driven by kinetic energy only. 
We also remark that we have, as usual, a (— ) branch, t < 0, with decreasing a and a future curvature 
singularity at t = and a (+) branch, t > 0, with increasing a and a past curvature singularity at t = 0. 
Our subsequent results will apply to both branches although, for the concrete discussion, we will mostly 
focus on the positive-time branch, where the universe expands. We find it convenient to use the scale 
factor a, rather than t, as the time parameter in the following. The remaining evolution equations can 
then be written in the form 

'A^" + rS^<A-^V^' = 0, (6.7) 

\gij4>''4>'' = 3, (6.8) 

where the prime denotes the derivative with respect to a. Hence, the scalar fields c/)^, viewed as 
functions of the scale factor a, move along geodesies in moduli space, with initial conditions subject to 
the constraint (|6.8D . 

Let us now apply these equations to the moduli space metric for the kink in a double-well potential, 
as computed in the previous section. To keep the formalism as simple as possible we will focus on the 
case of a small kink width, that is, e~^/U ^ 1. The moduli-space metric is then specified by Eqs. ( |5.9D , 
( IS.lOf ) and ( |5.8| ). Inserting this metric into Eq. ( |6.7D we find 

</." + eoe^-*Fz'^ = 0, (6.9) 

/?"-^eoe^-'^F(l + eV'^)^'' = 0, (6.10) 

+ (/3' - </>')^' + eV^^/3V - ^eV^i^'^ = 0, (6.11) 
while the constraint (16.81) turns into 



i<^'V |/?'' + ^eoe^-'^i^^'' = 3 . (6.12) 

The functions K = K{P, z) and L = L{P, z) are related to derivatives ol F = F{(5, z) and can be defined 
in terms of Jq, Eq. (^), as follows 

F{(i, z) = Jo (eVHl - z)) - Jo (-e^^^^^) , (6.13) 

, , {l-z)jUe^fi-\l-z))+zJ'(-ePfi-^z) , , 

m.z) = ^ ' ^ (6.14) 
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Figure 3: The function K and L which enter the effective equations of motion for the kink as a function 
of z for e^fi'^ = 10. 

The typical shape of F has been indicated in Fig. Fig. ^ shows the shape of K and L as a function 
of z. 

The equations of motion are generally quite complicated due to these functions. However, as the 
figures show F, K and L are non-trivial only in small regions around the boundaries with size set by 
He~^ (the width of the kink relative to the orbifold size) while they are relatively simple outside these 
critical regions. It is, therefore, useful to discuss the asymptotic form of the equations of motion away 
from the boundaries. First of all, for z G [0, 1] and away from the boundaries we have 



F ~ 1 



K ~ 



L ~ 



(6.16) 



Hence, for the kink being well inside the physical space the equations of motion (|6.S|)~( |6T2 ) greatly sim- 



plify and become, in fact, identical to the analogous equations derived from the M-theory action ( 2.13 ). 
On the other hand, for z < and away from the boundary we have 



F ~ 



L ~ -4 



(6.17) 



There are analogous results for z > \ but we will focus on the case z < for concreteness. Inserting 
these asymptotic expressions, we see that the equations (|6^), ( 6.10| ) and (|6.12|) for (p and [3 decouple 
from the z equation. They become, in fact, the equations for freely rolling radii and can be easily 
integrated to give 



4) = 3p^a + (1)0 , (3 = 3p/3a + (3o 
20 



(6.18) 



where (j)o and Po are arbitrary constants and the expansion powers and pf^ satisfy the constraint 

4 



4 + 3p| = 3 (6-19) 

which follows from ( |6.12| ). The evolution of the kink can now be studied in the background of these 
freely rolling radii. Inserting the above solutions for and /3 into the equation for z, Eq. ( 6.11| ), we find 



z" + 36z' + IfiQ^e^PP'^iQpisz + z')z' = , (6.20) 

where 



/.o = 5^ (6.21) 



is the width of the kink relative to the orbifold size initially at a = and 

S = P13-P<^- (6.22) 

Hence, for z < and away from the boundary the evolution of the kink is described by the single 
differential equation ( |6.2[1| ). 

7 Kink dynamics and kink-boundary collision 

We should now study the solutions to the system ( |6.9D -( |6yT^ ). Given that our main interest is in the 
collision of the kink with a boundary, ideally, we would like to find solutions with z S [0, 1] initially 
which evolve towards z ^ 0. Given the complexity of the equations, we cannot possibly hope to achieve 
this analytically. Later, we will address this problem numerically. However, some progress can be made 
analytically as long as z is away from the boundaries by using the approximate equations for z £ [0, 1] 
or z < discussed in the previous section. One may hope that finding such analytical solutions for the 
evolution up to shortly before and after the collision will lead to a correct qualitative picture of the 
collision process, roughly by gluing together these two types of solutions across the critical boundary 
region. As we will see in our numerical analysis, this is indeed the case. 

Let us start by looking at the case z £ [0, 1]. As discussed above, as long as z is not too close to 
one of the boundaries, the equations of motion reduce to the ones obtained from the M-theory effective 
action ( p. 13 ). Their solutions have been found in Ref. |2S] and are explicitly given by 



i_ 

(j) = 3p^^ia + 3{p^j-p^,i)ln(l + e'^^^'') + <Po, (7.1 

P = 3p;3,i« + 3(p^j-p^,,)ln(l + e-3'5»°)'^+/?o, (7.2 

d 

^ ~ 1 + e 
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35-^ + ^0 • (7.3) 



Asymptotically, for a ±00, these solutions approach freely rolling radii solutions for </) and /3 while z 
becomes constant. The early (late) rolling radii solution is characterised by the expansion powers p^^j 
and "pp^i {P(i,j and Ppj)- Both sets of expansion powers are subject to the constraint 

Pln + ^Pln = \ (7.4) 

where n = i, f and are related by the linear map 




Further, we have defined the quantity 

Si = Pi3,i - P<f>,i (7.6) 
which can be restricted, without loss of generality, to 



di > (— ) branch, 

5i < (+) branch. 



(7.7) 



We remark that 5f, the analogous quantity at late times, is given by 

Sf = PPJ - P<I>J = (7.8) 



as follows from the map ( [7.5[) . The remaining integration constants (pQ, /3o, zq and d are subject to the 
restriction 

^o-Po = In ( ^) . (7.9) 



Note that zq specifies the initial position of z which moves by a finite coordinate distance d to its final 
position zq + d. 

What is the relevance of these solutions in our context? First, we remind the reader that the above 
solutions play a double-role as exact solution to the M-theory effective action ( p.l3|) and approximate 
solutions to the kink effective theory if z £ [0, 1] and away from the boundaries. In their former role 
they present another indication that the effective M-theory action ( ^.13 ), as it stands, is not adequate 



to describe the collision process since the boundary values z = 0,1 are in no way singled out. In other 
words, z, as described by these solutions, passes through the boundary without being effected at all. For 
this reason, they will also be very useful for comparison with solutions to the kink evolution equations, 
to explicitly see the boundary effect in the latter. In their role as approximate solutions to the kink 
evolution equations for z £ [0, 1] they tell us that the collision can be arranged or avoided depending on 
a choice of initial conditions. Indeed, the initial position zq of the kink and the coordinate distance d by 
which it moves can be chosen arbitrarily. Hence, for the choice zq £ [0, 1] and zq + d G [0, 1] (and both 
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values away from the boundaries) the entire evolution of the kink is described by the solutions above and 
a collision with the boundary never occurs. There is, however, a caveat to this argument. While the kink 



becomes static asymptotically also the strong-coupling expansion parameter e necessarily diverges [23], 
as can be seen from the above solutions. Therefore, we eventually loose control of our approximation 
and the effective theory breaks down. Clearly, from the arguments so far, we cannot guarantee that 
the kink remains static when this happens. In this paper, we will not attempt to improve on this, 
for example by going back to the five-dimensional theory. Instead, we will be content with arranging 
a certain characteristic behaviour, such as the kink becoming static, to occur for some intermediate 
period of time before we loose control over the effective theory. 

Let us now analyse the evolution of the kink for z < and away from the boundary (the case 
z > 1 is similar, of course). In this case, the system is adequately described by the single approximate 
equation ( |6.20| ) for z while (j) 'ind P are decoupled and evolve according to one of the rolling radii 
solutions ( 16.181 ). Unfortunately, we did not succeed in integrating the z equation in general. However, 
we can find a number of partial solutions which, as we will see, provide a good indication of the various, 
qualitatively different types of z evolution. 

Let us consider the evolution of z in the background of a special rolling radii solution with a static 
orbifold, that is, 

pp = 0, p^ = (7.10) 



where the two possible values of Ptj, follow from Eq. ( |6.19|) . The equation ( 6.20 ) for z then simplifies to 



~" 3p^z' + 2fif^^z'^ = . (7.11) 



z 



The general solution to this equation can be readily found to be 



Zo + i: — In 

3p<i> 



1 + (e3P0« _ 1) 



(7.12) 



where zq and vq are integration constants specifying the initial position and velocity of z at a = 0, that 
is, zq = z{a = 0) and vq = z'{a = 0). Here, we are interested in solutions where zq is negative and as 
close to the boundary as is compatible with the validity of (|7.1l| ). In addition, we need vq < so z 
evolves into the region well-approximated by (|7.11 ). The parameter Vc is defined as 



3 

Vc = -P0/^o • (7.13) 

Let us discuss the properties of this solution for an expanding universe starting with the case Pif, = 
+ 2/^/3. It is easy to see from Eq. ( |7.12| ) that, independent of the initial velocity vq, z always diverges 
to — oo at some finite value of the scale factor a, in this case. For p^ = —2/^/3, however, the situation is 
somewhat more complicated and depends on the relation between \vo\ and \vc\- One has to distinguish 
the three cases 
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1) l^ol < |fc| : ^ converges exponentially to a constant 

2) \vo\ = \vc\ ■ z diverges to — oo as a ^ oo 

3) l^ol > \vc\ ■ z diverges to — oo at a finite value of a. 

Hence, we see that Vc plays the role of a critical velocity. As we will confirm later, these three cases 
already represent the three types of qualitatively different behaviour which can be observed for the full 
z-equation ( |6.20D or even the complete system (|6.9| )-( 6?T^ ). 

We should remark, though, that the second case l^ol = \vc\ while typical in that z diverges as a — > oo 
is not representative as far as the nature of the divergence is concerned. While its divergence is linear 
in a, the more characteristic case is an exponential divergence in a. The existence of such exponential 
divergences can be seen from the special solution 



to Eqs. ( 6.20 ). While this represents an exact solution for all values of and Pj3 we have to restrict 
signs to Pjs < and > so that z is negative and moves towards — oo. Within this range of and 
Pl3, however, the above solution shows an exponential divergence of z as a — > oo. 

After having identified the qualitatively different types of z evolution we can now ask more system- 
atically, based on the z equation ( |6.20 ), which type is realized for a given set of parameters and initial 



conditions. As can be seen from a rescaling of z in Eq. ( |6.2[1| ) the type of evolution cannot depend 
on the value of ^o- The only possible dependence is, therefore, on (recall that, for given ptf, is 



determined, up to a sign, from Eq. (6.19)) and the initial velocity vq = z'(a = 0). A relevant question 
in this context concerns the stability of the solution z = const which can be viewed as the limit of the 
exponentially converging case 1. Writing 

z = zo + C(a) , (7.15) 
where zq < 0, the linearised evolution equation for is, from Eq. ( |5.20 ), given by 



C" = -3 [5 + Afi^hoppe^P^'^] C' ■ (7.16) 

We conclude that the solution z = const can only be stable if 

Pf3 <0 and 6 = pp - p^ > . (7.17) 

It is only then that we expect the first case of convergent z to be realized. 

This can indeed by verified by a numerical integration of Eq. ( 6.20| ). Solutions with converging z 



exist if and only if the conditions (|7.17 ) are satisfied and, in addition, if the initial velocity \vo\ is smaller 



than a certain critical velocity Vc- A simple scaling argument shows that 

Vc = h{pp,p^)^io (7.18) 
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where /i is a function which, from the numerical results, turns out to be of and slowly varying. 



What happens outside the region (7.17)? If we leave this range by crossing = we find for small 
positive p/3 and \vq\ below the critical velocity that z still converges at first but then, in accordance with 
our analytic argument, develops an instability, which drives it to —oo at finite a. The intermediate 
stable phase gradually disappears as one increases p^. For > and \vo\ above the critical velocity 
one always finds divergence to — oo at finite a. Hence, for > we are always in the third case above. 
As we leave the region ( [7.17 ) crossing 6 = we find case 2 is realized below and case 3 above the critical 
velocity. However, as 5 becomes more negative, the critical velocity decreases rapidly until we are left 
with case 3 only. 

In summary, the converging case 1 is only found in the range ( [7.17 ) and for initial velocities smaller 
than a certain critical value while otherwise z always diverges to — oo typically according to case 3 at 
finite scale factor a. 



We can now try to combine the information we have gathered about the evolution of the system 
before and after the collision to set up criteria which will allow us to decide the outcome of a collision 
process. Let us consider a particular solution (7.1)~( |7^ for the evolution inside the interval z G [0, 1]. 
As we have already mentioned, the distance by which the kink moves is a free parameter so a collision 
may never occur. Then, this solution describes the full evolution of the system as far as it is accessible 
within the four-dimensional effective theory. On the other hand, if initial conditions are chosen so that 
a collision does occur, the particular solution ( |7.1[) ~ ([7.3D will determine the velocities z^^j, 0^^^ and P'^^i 
right before the collision. We can then, approximately, identify vq ~ z'^^^, 3p,f, ~ (p'^^i and 3p/3 ~ f3'^^i and 
apply the previous results for the evolution at z < 0. One concludes that only for a very low-impact 
collision with small z^^j and an orbifold size which, at collision, decreases less rapidly than the dilaton, 
that is P'^^i < and (p'^^i — fi'^^i > 0, does z converge to a constant. Otherwise z diverges to — oo and 
this can, in fact, be viewed as the generic case. 



Of course, the criteria above may be somewhat inaccurate since we have ignored the complicated 
structure of the evolution equations near the boundary. We have, therefore, numerically integrated the 
full system ( [6.9| )~( [6yT^ ) to test the above criteria for the outcome of a collision process. It turns out 
that, in broad terms, the picture remains qualitatively the same. 

Starting with z near zero inside the [0, 1] interval, we went around the ellipse ((/'[.oi)^ -|-3(/3^qj)^ ~ 12. 
Note that in this case the exact identity cannot be observed since the constraint equation Eq. ( [6.12 ) 
includes an extra term proportional to (^^oi)^- Nevertheless the correction is always small since we set 
'/'col ~ Pcoi to a large negative value. This makes the initial value for e very small and allows us, for 
the cases where e grows, to follow the evolution for longer times until e ~ 1 and the four-dimensional 
effective theory breaks down. We also chose a large initial Pcoi so that e~^n remains as small as possible 
during the evolution, for the cases with /J^^j < 0. In all cases we set eo = 1 and fi = 0.2. 
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Figure 4: Position modulus z for the kink (solid line) and M-theory three-hrane (dashed line) as a 
function of the scale factor a. The initial conditions have been chosen as z^oi = 0.027, z'^^^ = —0.12, 
/?coi = 2.0, = -0.72, ,/.,oi = 16.14, cP'^^^ = -3.23 



For each of these sets of initial conditions we then varied z^^j from zero upwards and looked for 
changes in the large time behaviour of z. The numerical results were obtained by evolving Eqs. (|6.S|)- 
( |6.1lD using a fourth-order fixed step Runge-Kutta method. The accuracy of the method was checked 



by confirming that the constraint equation Eq. ( 6.12| ) was satisfied throughout the evolution. The 
individual terms on the left hand side of Eq. ( 6.12| ) should sum to 3, and typically after 2000 time- 
steps of size 0.01 the deviation from this value was smaller than 0.01%. In the worst cases, where the 
equations of motion are no longer valid because one of the assumptions have broken down, the sum 
never gets above 0.2% 

In Fig. ^ we have an example of the first type of behaviour, for a small negative value of (3'^^^ After 
crossing the boundary at z = the kink relaxes to a stable constant solution. For early times this 
solution matches the one obtained from the M-theory effective action for the same initial conditions. 
Nevertheless, as soon as the kink approaches the boundary the two start differing, converging to different 
asymptotic values. 

For a slightly higher value of the initial velocity the difference is even more striking, as shown in 
Fig. |5[ In this case z diverges in finite time, indicating that we are above the critical velocity. This third 
case turns out to be the most common, as already observed in the simplified system. Only for P'^^y < 
and /3^qj — cp'^^^ > and z^^j below the critical velocity does the system avoid singular behaviour. 
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Figure 5: Same as in Fig. but with z'^^^ = —0.14. 



In Fig. H we have an example for a solution corresponding to case 2. Here both P' , and (p' , are 



'^^i < 0, z does not 



negative and we are below the critical velocity. As a consequence of 5 ~ (3'^^^ 
relax to a constant but its magnitude increases exponentially instead. In this case the solution has to 
be taken with care, since e~^/i quickly becomes large in the exponential regime and the equations of 
motion stop providing a reliable approximation. 

Finally we have checked that once we go above the critical velocity, z always diverges for finite a. It 
is well known that in cp'^ theory when a kink anti-kink collision takes place, above a certain limit velocity, 
they reflect and bounce back (for lower velocities they can either reflect or form a bound state). 
This behaviour is a consequence of a resonance effect between the kink pair and higher field modes, 
so we should not be surprised not to observe it in the context of our four-dimensional effective action. 
This does not yet exclude the possibility of a bounce in a high-velocity regime which is accessible only 



in the context of the full five dimensional theory, a question which is currently under investigation |28]. 



What do these results imply in terms of the five-dimensional defect model ( |3.1D ? As we have seen, 
if z starts its evolution within the interval [0, 1] and subsequently collides with a boundary (at z = 0) it 
is generically driven to — oo very rapidly. It should be stressed that the z kinetic energy remains finite 
at this singularity. Nevertheless, we do expect the effective four-dimensional theory to break down 
eventually, as z — oo. This is because some of the higher-order term we have neglected are likely 
to grow with z, in a way similar to the linear z term in Eq. ( |6.20| ). However, at least for sufficiently 
small expansion parameters e and ne~^ the four-dimensional theory will be valid some way into the 
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Figure 6: Position modulus z for the kink (solid line) and M-theory three-brane (dashed line) as a 
function of the scale factor a. The initial conditions have been chosen as Zcoi = 0.027, z'^^^ = —0.060, 
Pcoi = 2.0, = -1.77, (^eoi = 14.75, </.^„i = -1.61 

singularity. Hence, wc can conclude that a five-dimensional kink, interpolating between the vacua 
X = nv and X = ('^ + 1)"^ which collides with the boundary at z = effectively disappears and leaves 
the field x the vacuum state x = + 1)'*^ (and an analogous statement holds for collision with the 
boundary at z = 1). Prom the M-theory perspective, such a process corresponds to a transition 

{PuP2,P3) = (n, -{n + 1), 1) {Pup2,p3) = (n + 1, -(n + 1), 0) (7.19) 

between two different sets of charges and, hence, topologically different compactifications. 

8 Conclusion and outlook 

In this paper, we have presented a toy defect model for five-dimensional heterotic brane-world theories, 
where three-branes are modelled by kink solutions of a bulk scalar field x- We have shown that the 
vacuum states of this defect model correspond to a class of topologically distinct M-theory models 
characterised by the charges Pi and P2 on the boundaries and the three-brane charge P3. Specifically, 
we have seen that a state where x equals one of the minima x = Xn = '^^ of the potential, where n G Z, 
corresponds to a state with charges (/3i, /32, /^a) = {n, —n, 0), that is, an M-theory model without three- 
branes. If, on the other hand, x represents a kink solution interpolating between the minima x = ^^'f 
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and X = ('^ + 1)^ the associated M-theory charges are {Pi, P2, Ps) = {n, — (n + 1), 1) corresponding to a 
model with a single-charged three-brane. 

We have computed the effective four-dimensional action associated to a kink solution and have 
studied the time-evolution of a kink in this context. Our results show that, generically, a collision of 
the kink with a boundary will lead to a transition between the two types of vacua mentioned above. In 
other words, the kink will disappear after collision which corresponds to a transition between a state 
with a single-charged three brane and a state without a three-brane. 

There are several interesting directions which may be pursued on the basis of these results. Clearly, 
our original M-theory model as well as the associated defect model are rather simple and a number 
of possible extension and modifications come to mind. First of all, we may try to modify our defect 
model by including more than one additional bulk scalar field, in particular to allow for exact BPS 
multi-kink solutions. One may ask whether the defect model can be embedded into a five-dimensional 

= 1 supergravity theory as is the case for the original M-theory model. Further, there are a number 
of generalisations of five-dimensional heterotic M-theory, such as including a more general set of moduli 
fields IQ, which one may try to implement into the defect model. For example, including the general set 
of Kahler moduli would allow one to study topological transitions of the underlying Calabi-Yau space 
through flops, in addition to the types of topology change considered in this paper. 

Perhaps the most interesting direction is to study the evolution of more complicated configurations 
of our defect model (p.l|). For example, one could envisage evolving the field x from some initial (say 



thermal) distribution to see which type of brane-network develops at late time P8|| . In particular, one 
would like to answer the important question whether the system can evolve from a brane-gas to a brane- 
world state. If this is indeed what happens such an approach will lead to predictions for the late-time 
brane-world that has evolved, given a certain class of plausible initial states. Concretely, within the 
context of the simple model presented in this paper, we may expect predictions for the charges Pi in 
this case. As we have discussed, the values of these charges are correlated with important properties of 
the theory such as the type of gauge group. Optimistically, we may therefore hope that our approach 
leads to prediction for such important low-energy data, at least within a restricted class of associated 
M-theory compactifications. 
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